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Abstract 

We apply our systematic NLO small x resummation of singlet splitting functions to 
the scaling violations of structure functions and compare the results with data. We develop 
various theoretical tools which are needed in order to relate resummed parton distributions 
to measurable structure functions, and we present results from a variety of fits to HERA 
data for the structure functions F2 and Fl using the resummation. The behaviour of 
the singlet splitting functions at small x and fixed Q 2 is effectively parametrized as x~ x . 
We find that, for A small or negative, the resummed description of scaling violations may 
be phenomenologically as good as or even better than the standard next-to-leading order 
treatment. However, the best fit gluon density and value of a s can be significantly modified 
by the resummation. 
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1. Introduction 

A complete understanding of scaling violations at small x in deep inelastic structure 
functions within the framework of perturbative QCD remains elusive. Whereas techniques 
for the inclusion of small x contributions to leading twist evolution equations have been 
known for some time [1,2], their precise implications have been rather unclear. On the 
phenomenological side, while the corrections should be sizeable, the next-to-leading order 
(NLO) description of HERA structure function data is so successful [3] that there seems 
to be little room for improvement [4]. On the theoretical side, the evaluation of small x 
corrections to the singlet splitting function using the BFKL theory [5,6] appears to fail: the 
recent calculation [7-10] of the NLO contribution xi to the kernel x = a^Xo + c^Xi--- shows 
that the expansion is very badly behaved, as the non leading term completely changes the 
qualitative features of the leading order. However, recent studies have led to considerable 
progress: in fact, it turns out that by assuming the simultaneous validity of leading twist 
Q 2 evolution and the BFKL equation, most of the theoretical problems encountered in this 
approach can be overcome. 

At large Q 2 and not too small but fixed x the QCD evolution equations for parton 
densities [11] provide the basic framework for the description of scaling violations. The 
complete splitting functions have been computed in perturbation theory at order a s (LO 
approximation) and a 2 (NLO) [12]. For the first few moments the anomalous dimensions 
at order a 3 are also known [13]. However, at sufficiently small x the computation of the 
splitting functions based on only the first few terms of the expansion in powers of a s cannot 
in general be a good approximation. Even assuming that a leading twist description of 
scaling violations is still valid in some range of small x, as soon as x is small enough that 
a s £ ~ 1, with £ = \ogl/x, all terms of order cn s (a; s £) n and a 2 (a s £,) n which are present [6] 
in the splitting functions must be considered in order to achieve an accuracy up to terms 
of order a^. 

In most of the kinematic region of HERA [3] the condition a s £ ~ 1 is indeed 
true. Since a s £ < a s \ogs/Q 2 , and s m 10 5 GeV 2 , a s (m 2 ) 0.119, we see that at 
Q 2 = 3, 10, 10 2 , 10 3 GeV 2 , a s £ can be as large as 4.3, 3.0, 1.2, 0.6, respectively. So one 
would expect many terms of the series to be important and, consequently, to see in the 
data indications of significant corrections to the approximation [14,15] in which only terms 
up to order a 2 are kept. In reality this appears not to be the case: the scaling violations 
in the data are in excellent agreement with the predictions of the leading twist evolution 
equations in the NLO approximation [3]. Of course it may be that some corrections exist 
but that they are hidden by a modification of the fitted gluon, which is the dominant parton 
density at small x, or by a change of the measured value of a s . However, a straightforward 
inclusion of small x corrections completely spoils [4] the agreement between the data and 
the NLO calculation. 

From a theoretical viewpoint, the inclusion of contributions to perturbative QCD evo- 
lution to next-to-leading log x presents two main problems. The first is that the qualitative 
features of the BFKL evolution kernel at leading- and next-to-leading log x are completely 
different. In particular, leading log x evolution leads to a parton distribution which rises as 
a power when x — > at fixed Q 2 [5] , while at the next-to-leading level the structure func- 
tion displays an unphysical oscillatory behaviour [16]. This is because at next-to-leading 
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order the function x(M) has a qualitatively different shape as a function of M (the vari- 
able which is conjugate to logQ 2 upon Mellin transformation): in particular near M = 
Xo ~ while xi ~ — 1/M 2 . Furthermore, even if this problem is side-stepped by 

treating the subleading correction near M = \ as a perturbation, there is a new difficulty, 
namely, that the ensuing correction to the asymptotic behaviour is very large [17,18]. 

The first problem is related to the presence of unresummed large logs of Q 2 in the 
leading In a; resummation [19-21]. It was recently shown by us [22] that it is possible to 
implement a reorganization of the perturbative expansion such that large logs of Q 2 and - 
are both resummed simultaneously (double- leading expansion [1]). This can be achieved 
by exploiting the fact that when Q 2 and ^ are both large, the perturbative evolution 
admits a dual description [22,19], either in terms of an evolution equation in t = lnQ 2 /A 2 
(the usual leading twist evolution equation [11]) or in terms of an equation in £ = lnl/a; 
(the BFKL equation [5]): leading and next-to- leading large logs can be determined from 
the knowledge of the fixed order kernels in both equations, and then used to construct 
double-leading kernels appropriate to either equation. 

It should be clear that our function %(M) which corresponds by duality to the leading 
twist anomalous dimension order by order in perturbation theory is not precisely the same 
as the Mellin transform of the kernel of the £ evolution equation: it cannot be since 
beyond LLx the latter is in general a differential operator. However the two are closely 
related. At LLx they are equal. If a s does not run, they are simply related to all orders 
through a change of scale [7]. When the coupling runs, duality still holds, however the 
relation becomes more complicated: it may be specified order by order at NLLx, NNLLx, 
etc. [23], but to all orders becomes ambiguous. What we do in the following is to extract 
unambiguous information from the NLLx BFKL kernel in order to learn about %(M) and 
thus by duality about the small x singularities of the leading twist anomalous dimension. 
Similarly, we use the NLO anomalous dimension to infer, by duality, information about the 
small M singularities of %(M). Combining these two independent sources of information, 
we are able to construct the double-leading expansion in which all large leading and next-to- 
leading logarithms are correctly resummed at leading twist into an anomalous dimension: 
logarithms of x at NLLx, logarithms of Q 2 at NLLQ through the usual renormalization 
group improvement. 

Since logarithms of Q 2 are resummed, the double-leading expansion turns out to 
be free of qualitative instabilities, thereby providing a complete solution [22] of the first 
problem mentioned above. The physical origin of this stability is the powerful constraint 
of momentum conservation, which is automatically taken into account in the double- 
leading expansion. By duality, this constraint fixes the all-order value of the kernel %(M) 
at M = 0, thereby stabilizing the expansion of the kernel in the neighbourhood of this 
point. At this stage one already understands why the corrections at small x are not 
catastrophic. It is important to notice that the possibility of matching evolution in InQ 2 
and evolution in lnl/a; [22] shows that there is no breakdown of factorization at small x 
and that the assumption of the dominance of leading twist terms is tenable, in agreement 
with the conclusions of refs. [24], based on model calculations. Our solution of the first 
problem is similar to that of ref. [25] which was obtained from a resummation of the 
BFKL kernel, but we think that our approach, which also includes all the information 
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contained in the conventional NLO leading twist evolution equation, is more general and 
avoids model dependent assumptions. Alternative approaches to this problem have also 
been suggested [26]. 

The second difficulty mentioned above arises due to the fact that the asymptotic be- 
haviour at small x and fixed Q 2 (Regge limit) of the perturbative evolution is subject 
to large corrections because nonleading terms of order a s (Q 2 )(a s (Q 2 )^) n are not asymp- 
totically small for x — > at fixed Q 2 . The effect of these terms is large and makes the 
small x expansion unstable. However it was recognized [20,27,23] that this instability is 
due to the fact that higher order contributions to the structure function must change the 
asymptotic small x behaviour from x~ x ° to x~ x = x~ x "e Ax ^ ~ x~ x "[l + AA£ + ....], 
with Ao = ol s xo(\) and AA = a 2 xi{\) + ■ ■ ■■ This problem can be cured by resumming 
the contribution to the asymptotic behaviour, i.e. effectively treating the whole of A as 
a leading order term. Indeed, in ref. [23] it was proven that this can be done order by 
order, thereby removing the uncontrolled growth of the subleading small x contributions 
mentioned previously. This procedure requires knowledge of the asymptotic behaviour in 
the Regge limit to all orders. However the true behaviour is quite possibly nonperturbative 
and, in any case, certainly cannot be reliably evaluated by just a few perturbative terms 
(the NLO term is larger than the LO one for realistic values of a s (Q 2 )), even though it can 
of course be computed within specific models [24-26]. We prefer to take a more general 
approach, and thus, at fixed Q 2 , we treat A as a parameter to be fixed from the data. 

For a given value of A one obtains for the double-leading resummed splitting functions 
a well behaved perturbative expansion which resums all leading logs of x and Q 2 , and 
should thus give an accurate description of scaling violations in a wide region of x and 
Q 2 . There is still a residual ambiguity, due to the fact that, in order to avoid double 
counting, in the double-leading expansion one must subtract order by order terms which 
contain leading logs of both x and Q 2 , and would thus otherwise be included twice. Since 
the small x part of the double-leading expansion has been stabilized by the resummation, 
there is an ambiguity in the treatment of these double-counting terms, in the sense that 
one can choose whether or not to include them in the resummation. 

We will see that that for reasonably small values of A (A ~ 0.2) the resummed struc- 
ture functions are quite close to the standard two-loop results, especially if the remaining 
ambiguity is solved in such a way that resummed splitting functions are as close as possible 
to the two-loop ones. Thus, the success of the two-loop description of structure function 
data from HERA can be understood. However, even in this case, we still find that the 
resummation can have a significant impact on the extraction of parton distributions and 
of the value of a s from the data. Thus one important practical conclusion of our work is 
that the gluon and the strong coupling extracted from the NLO fit to HERA data at small 
x are to some extent biased. The precise size and direction of this bias depends on the 
particular resummation adopted and on the value of A. This introduces an extra source 
of uncertainty in parton distributions and on the value of a s extracted from scaling viola- 
tions, which must be taken into account for a correct estimate of the associated theoretical 
error. 

It is the purpose of this paper to discuss the phenomenology of structure functions 
at small x on the basis of the treatment of splitting functions developed in ref. [22] . Our 



3 



aim is to provide a comprehensive self-contained treatment, which on the one hand gives 
all the technical details which are needed for a practical implementation of the resumma- 
tion procedure, and, on the other hand, discusses several theoretical issues raised by our 
resummation method. As discussed above, our resummation consists of two steps: first, 
one constructs the double-leading expansion of anomalous dimensions and their associated 
splitting functions, and then, the small x asymptotics is resummed thereby introducing a 
dependence on the parameter A. Likewise, the resummation of structure function proceeds 
in two corresponding steps. First (sec. 2), we construct the double-leading expansion for 
the physically observable structure functions F<i and Fl- Then (sec. 3), we resum the 
ensuing expression of the structure function by extracting the asymptotic behaviour in the 
Regge limit. The choice of factorization scheme and the theoretical ambiguities which it 
induces in the resummation procedure and running of the coupling, are discussed in sec. 4. 
Readers who are only interested in the construction of the resummed structure functions 
in the DIS scheme can skip sec. 4, and turn to sec. 5, where we fit recent HERA data on 
the F 2 and Fl structure functions within our framework, and show that the resummation 
somewhat improves the quality of the fit and its stability and has a significant impact on 
the determination of the gluon density and strong coupling constant. 



2. Double— leading expansion 

In our previous work [22] we have derived a resummed expression for the largest 
eigenvalue of the singlet anomalous dimension matrix and its associated splitting function 
in the double-leading expansion. In this section, we construct the double- leading expansion 
of the full matrix of anomalous dimensions and coefficient functions, and construct the 
resummed structure functions F 2 and Fl in the DIS scheme, which we will later use for 
the fits to data. Other schemes will be discussed in sec. 4. We will henceforth adopt the 
notation, conventions and terminology of ref. [22]. Even though we will recall all results 
which are needed for the resummation of the structure function, we refer to ref. [22] for a 
more complete derivation of the resummation at the level of the largest eigenvector. 



2.1. Double leading anomalous dimensions 

First, we wish to discuss the full two by two matrix of anomalous dimensions. Recall 
that [6] only the gluon entries of the anomalous dimension matrix contain leading log x 
(LLx) singularities, i.e. 

-( ^ ou 

7LLx I gq I a s \ „,qg {ce s \ ? \ Z - 1 ) 



where the nonvanishing entries satisfy the color-charge relation 



7f = 7f , (2-2) 



where the color factors are Cp = % — and Ca = n c . This implies that at this order the 
eigenvalues of the anomalous dimension matrix are 7+ = 7f s and 7" = 0. Because the 
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small eigenvalue vanishes at the LLx level, it is possible to set it to zero to all orders by 
choice of factorization scheme [28,19]. 

At the next-to-leading log x (NLLx) level the eigenvalues are given by [1] 



ry9Q 

7^ LLx = 7f + a, ( 7 f| + j^lfs) + 0{a 2 s ) 

Is 

ry9Q 

7N LLx = « s (7^-^7|f)+0(^). 



(2.3) 



while the corresponding eigenvectors are 



'19 



Q+ = a s ^§G+ +0(a 2 s ), 
Q- = - 2 j m G- + 0(a s ). 



(2.4) 



It is easy to see that a sufficient condition for singular contributions to 7 to vanish at 
this order is that a color-charge relation also holds in the quark sector, i.e. that 



7. 



qq _ Cf_( qg 
— c A ^ ' ss 



eg), 



(2.5) 



where e\ is a scheme dependent constant. This condition is respected in both the MS and 
DIS schemes [29], and we will henceforth only consider factorization schemes in which it 
is satisfied. 

Since in any such scheme the small eigenvalue 7" is regular at N = 0, and therefore 
unaffected by the summation of leading log x singularities, it is possible [17,18] using 
eq. (2.3) to reconstruct the full matrix of anomalous dimensions to next-to-leading log x 
from the knowledge of the NLLx large eigenvalue [7], and of the NLLx quark anomalous 
dimension [29]. However here we are interested in the double-leading expansion of the 
evolution equations. In this expansion, the anomalous dimension is constructed by adding 
the usual one and two loop contributions on top of the leading and subleading singularities, 
and subtracting the double counting [22]: 



7 + L (iV,a s ) = [a s7o + (A0 + 7 S + (ft) - 



+ a s 



a s7l + (A0+7 s + s (f)-« 



s AT 



,+ 



+ 



(2.6) 



where 7q~ and 7^ are the one and two loop contributions to the largest eigenvalue, 7+ 
and 7+, are the leading and subleading singularities, and the double-counting subtractions 
in the DIS and MS schemes are given by ef = nfn c (5 + 13/(2n^))/(187r 2 ) and = 
— (lf n c + n f)/(^ 1Tn c)- The small eigenvalue 7" instead simply coincides with its two loop 
form: 



7 DL (iV,a s ) = a s7o (N) + a 2 s Yi W 
Likewise, we can construct a double-leading anomalous dimension matrix as 



7 g L (iV,« s )= « s7 f(Af) + 7 ^(ft)-d.c. +a s a sl l J (N)+^ s {^)-d.c. 



+ • 



(2.7) 



(2.8) 
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where i,j = q,g and double counting term are subtracted as in (2.6). It is easy to prove 
that the eigenvalues of this double-leading anomalous dimension matrix are given order by 
order in perturbation theory by the double-leading sum of the eigenvalues a s ^-\-a 2 s ^-\-. . . 
of the large x evolution matrices « s 7q J + a^- 7 + . . . and eigenvalues 7^ + a s 7^, + ... of the 
small x evolution matrices 7J 7 +a s 7*^ + . . ., after subtracting the double counting terms in 
each case. However because the eigenvalues depend nonlinearly on matrix elements, this 
linear relation between the double-leading eigenvalues and the eigenvalues of the small 
x and large x evolution matrices only holds order by order up to subleading corrections. 
Therefore, it is in practice more convenient to use as primary quantities the eigenvalues 
(2. 6), (2. 7) and two matrix elements, and determine the other two matrix elements from 
them. 

We thus choose to adopt the eigenvalues and quark sector matrix elements 7" and 
7" as primary quantities. This choice is motivated by the fact that the large eigenvector 
is the primary quantity in the small x evolution equation (BFKL equation) whose kernel 
determines the small x contributions to the anomalous dimension 7+ by duality [22] ; while 
the quark sector small x anomalous dimensions are directly determined through all-order 
small x factorization (kx factorization) of the deep-inelastic structure functions F 2 and 
F L [29]. In double leading expansion 

1 q D i(N,a s )=a s ^ 9 (N) + a s 
^i(N,a s )=a sl q \N)+a s 

where (in DIS scheme) = n//67r, e\ = 13n c n//367r 2 , and in the second expression we 
have used the color-charge relation (2.5). Given the eigenvalues and quark entries we can 
then determine 7 s9 from the trace condition, and 7" from the determinant: 

7 + + -y- = Y 9 + 7", 7 + 7~ = 7 33 7 <?<? - 7 g3 7 ffg . (2.10) 

Note however that while to determine r ) gg at NLLx it is sufficient to know the eigenvalues 
and r ) qq at NLLx, to fix 7" from the determinant condition would require 7" and r ) qg 
at NNLLx. It follows that the value of r ) 9q at NLLx is of no consequence for a NLLx 
calculation. 

As already discussed in ref. [22], at next-to-leading order in the double leading ex- 
pansion, momentum conservation will be violated by next-to-next-to leading terms, and 
can thus be restored by adding a subleading correction. Indeed, momentum conservation 
implies that, at N = 1, a s ) + 7^(1, a s ) = and 7 9S (1, a s ) + 7 S9 (1, a s ) = 0, so 

7 +(l, a a ) = 0, 7 "(1, a a ) = 7 M (1, a s ) - 7 93 (1, <*,). (2.11) 

While in DIS or MS the one and two loop contributions to the large eigenvalue 7^ L eq. (2.6) 
already vanish at N = 1, the singular contributions 7 S and 7 SS give a non- vanishing 
contribution of 0(a^). Since this violation is sub-subleading we are free to remove it by 
subtraction, i.e. by setting 

7dl(W. <x 8 ) | mom = 7 + L (W, a 8 ) - t£ l (1, <**)■ (2-12) 

6 



+ 



*sll q {N) + % 



Iss \ N ) 



a 



: N ^OJ 



(2.9) 



+ 



Likewise, the condition on the small eigenvalue 7" is automatically respected by the one 
and two loop contributions to 7^ an d 7dl ec l- (2-9), but is violated by the singular terms 
7|f : recalling that 7" is free of N = singularities, one sees that the condition on 7" 
eq. (2.11) would require that, when N = 1, 7ff (a s ) = a s e\ + eg. So the small eigenvalue 
condition is also violated by terms of O(a^), and can be restored by a sub-sub leading 
subtraction: 

7dl(- /V > a s )lmom = 7dl(- /V > ««) - «»(7«(a 4 ) ~ a s e i ~ e o) 

7 gL(JV, a s )| mom = lH(N, a.) - g« s ( 7 |f(« s ) - « s e? - eg), ^ 

consistent with the color-charge relation (2.5). 

The anomalous dimension matrix to NLO in the double-leading expansion is fully 
determined by eq. (2.10), using the double-leading eigenvalues eqs. (2.12), (2.7) and the 
quark matrix elements eq. (2.13). All the quantities which are needed for a NLO compu- 
tation are known. Specifically, the NLO singularities of the large eigenvector 7+ can be 
determined using the NLO duality relations 

TV 

Xo(7+(f )) = -, (2-14) 

OLa 



its(^) = - ]y_ s + )i[[ , (2.i5) 



Xi(7+(f)) 

in terms of the well-known BFKL [5] kernel 

X0 (M) = [V>(M) + V(l - M) - 2^(1)] (2.16) 

7T 

and the NLO kernel %i, which was determined in refs.[7-10] in a scheme which is closely 
related to MS, and in the MS scheme is given by [23] 

Xi (M) = ^n 2 J(M) + ^/3 n c (2V/(l) - ip'(M) - V>'(1 - M)) + i^Xo(M) 2 , (2.17) 

where the function 5 is defined in the first of ref. [7], and (3q = ^-n c — fro/. As we will 
show below, the same expression also holds in the DIS scheme. The quark singular terms 
list lis have been computed in the MS and DIS scheme in ref. [29]: in DIS 

*sl s q s g m = h 2 (ls(^))R(ls(^)l (2.18) 
where h 2 (M) is a process dependent contribution (given by eq. (5.20) of ref. [29]), while 
R(M) is the process independent gluon normalization factor (given in ref. [29] as the 
solution (3.17) of the differential equation (B.18)). 

In practice the anomalous dimensions are evaluated by first computing the coefficients 
in their expansions in powers of jf: 



00 



^ s ( N ) ~~ 41n2 ^2 ttn (/v) ' 
n=l 

00 00 

w*) = (1 + E * (t r ) - a (1 + e "'n »r ) (s-w) 



n- 

00 



^19, <± 



(t) = S(i + E c »(?f)") 

n=l 
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where a s = 41n2n c a s /7r, and the coefficients are thus normalised so that each series has 
radius of convergence unity. The corresponding series for the splitting functions are uni- 
formly convergent for all x > [1]: enough coefficients for sufficiently precise calculations 
in the HERA kinematic region are given in tables 1 and 2. 

To determine the evolution of the singlet quark and gluon by solution of the singlet 
evolution equations an explicit determination of the complete two by two matrix of anoma- 
lous dimensions is in fact not necessary. Rather, we can conveniently construct the solution 
from the eigenvalues and quark-sector entries by means of a projector formalism [30] . The 
anomalous dimension matrix is decomposed as 

7 = M+7+ + M_7" (2.20) 

where 7 ± are the eigenvalues of the matrix 7, and the projectors satisfy 

M + +M_= 1 ; M±M± = M±; M + M_ = 0. (2.21) 

Explicitly, the projectors are given in terms of the eigenvalues and quark-sector entries of 
the anomalous dimension matrix by 

M - 1 ( ^11 ~ 7" ll9 \ . 

^ / 7+ " 7M ' <**> 

M_ = - ( 7 " T« _ ] 

7+ - 7" V - X 7g 9 - 7 ; ' 

where X = (7+ - 7<?<?)(7 ?9 ~ l~)/l qg - The evolution equation j-J = ~/(N,a s (t))f for 
the vector f(N,t) = (^^j) is then immediately solved in terms of the path-ordered 
exponential 



f(N, t) = Pexp f dt' [M+(N, a s (t>))<y+(N, a s (t')) + M_(iV, a s (t / ))7"(^, a s (t'))] f(N, t ). 

J to 

(2.23) 

Introducing the double-leading expansions (2. 6), (2. 7), (2. 9) of 7" and 7 9S into the 
eq. (2.22) we thus obtain a double-leading expansion of the projector itself, which in turn, 
once substituted in eq. (2.23) gives the standard next-to-leading solution [30], but now in 
the double-leading expansion. Indeed, for numerical computations the solution eq. (2.23) 
can be used directly. This completes the construction of the double-leading approximation 
to the evolution equations for parton distributions, since the nonsinglet and valence quark 
distributions are free of small x logs and can be treated at two loops in the usual way. 

2.2. Double leading structure functions 

We can thus proceed to the determination of the structure functions. In the DIS 
scheme [31], the structure function F2 simply coincides with the quark distribution: 

F 2 (x,t) = (e 2 )2n f Q(x,t) + F^ s (x,t), (2.24) 
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where (e 2 ) = ^^e, 2 , Q(x,t) = £ E"ii (ft(M) + 9i(M)) and F 2 NS (x,t) is the nons- 
inglet component. However, the identification of F2 with the quark distribution still leaves 
some freedom in the definition of the gluon distribution. Introducing the general scheme 
change matrix Uij(N, a s ), with i,j = q,g, such that (in matrix notation) 

f'(N,a s ) = U(N,a s )f(N,a s ), (2.25) 

, qDIS . , QMS , nT q 

then if / = ( ^ DIS J and / = ( g M s ) > the condition F 2 = Q only determines the quark- 
sector matrix elements U qq = C\ (N, a s ) and U qg = Cf (iV, a s ), where C l 2 are the quark and 
gluon F 2 coefficient function in the MS scheme, leaving the gluon sector matrix elements 
U gq and U gg undetermined. This freedom can be fixed by assuming [32] the validity 
for all moments of the relations (2.11) which momentum conservation imposes on second 
moments, i.e. Ug g {N,a s ) = 1 — Cf (iV, a s ), U gq (N,a s ) = 1 — Cf(iV, a s ). Because the 
MS coefficient functions are free of leading singularities, with this choice the singular 
contributions to the eigenvalues of the anomalous dimension matrix up to NLLx are the 
same in the DIS scheme as in the MS scheme. Furthermore, the Fl coefficient functions, 
which start at next-to-leading lnx, are also the same in the MS and DIS scheme. Since 
the eigenvalues and quark entries of the anomalous dimension matrix fully determine the 
structure function, the scheme is fully determined by this somewhat weaker assumption, 
which we adopt as a definition of the DIS scheme in the double leading expansion. We 
will come back to a fuller discussion of the relation between double-leading MS and DIS 
schemes in sec. 4. 

We can thus determine easily the double-leading expansion of the structure functions 
in terms of the double-leading parton distributions: F2 is just given by eq. (2.24), while 
Fl is constructed from Q and G using double-leading coefficient functions. These are in 
turn constructed as the sum of two loops and leading singular contributions, minus the 
double counting: 

CI, DL (iV, a s ) = a s [CI ,{N) + ss (%) - e\] , (2.26) 

where i = q, g. The two loop coefficients i(N) were computed in ref. [31], and the 
singular terms ss (^f) in ref. [29]: 

*sCl ss (ft) = h L M^))R(ls(^)), (2.27) 

where the process dependent piece /il(M) is given by eq. (5.14) in ref. [29], and C q L is given 
by the color-charge relation ss = ^(Cl ss — e 9 L ). The double-counting subtractions 
are e 9 L = e q L = 0, and the coefficient functions are again best evaluated through their 
series expansion 

00 

CfU(t) = ^(l + £^(f) n ), (2.28) 

n=l 

where the coefficients dk are listed in table 2. 
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3. Resummation of the structure function 

A resummation of the expansion of the splitting function at small x is required [22,23] 
in order to obtain a stable expansion in powers of a s . The resummation affects the small 
N behaviour of the anomalous dimension, and specifically the expansion of the anomalous 
dimensions in leading, subleading,. . .singularities: [1,22]: 7 = [y a + a s 7 ss + ...]. In this 
section we will discuss how this resummation affects the determination of the structure 
functions. The resummed version of the small x expansion will then finally be combined 
with the usual loop expansion in order to construct a resummed double-leading expansion. 

3.1. Resummed anomalous dimensions 

The resummation is based on treating the asymptotic small x behaviour of the splitting 
functions as effectively leading order. This means that a constant is subtracted from the 
contribution to the BFKL kernel at each perturbative order, and added to the leading 
order: 

X (M, a s ) = a sX o(M) + a 2 sX i(M) + ... (3.1) 
= a sX o(M) + a 2 sX i(M) + ..., (3.2) 

where 

a sXo {M,a s ) = a sX o(M) + AX(a s ), Xi {M) = Xi (M) - a (3.3) 
for i = 1,2,..., and thus 

00 

AA(a s ) = ]T< +1 Cn . (3.4) 

n=l 

The constants q are uniquely fixed [23] order by order by the requirement that the asso- 
ciated splitting functions P+, -P+,..., define a stable expansion, in the sense that 

lim P ss (x, a s )/P s (x, a s ) = f{a s ), (3.5) 

where f(a s ) does not depend on x, and thus in particular does not grow when x decreases. 
At next-to-leading order the constant is simply equal to the value of the subleading cor- 
rection to the BFKL kernel evaluated at the leading-order minimum: 

Ci = Xi(i)- (3-6) 

More generally, if X (M) has a minimum then Ci are such that the minimum of X o(M) 
coincides with the minimum of X . For convenience we define Co = Xo(^) 5 so that if 
X (M, a s ) has a minimum, 

00 

A(oj s ) =a s ^2 a ™°n ( 3 - 7 ) 

n=0 

is its value at that minimum. 

The sum in the definitions of A and AA is to be understood as a symbolic indication 
that even though AA can be formally expanded order by order in perturbation theory, it 
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is its all-order value which determines the behaviour of the splitting function in the Regge 
limit, and is thus relevant for phenomenology at small x. It is important to notice that 
this is an inevitable consequence of the assumption that the coupling runs with Q 2 , and 
it is thus common to any perturbative computation based on this assumption, such as 
those of ref. [25]. Of course, a fixed order computation may turn out to provide a good 
approximation to the all-order value of A. However, the known terms suggest at best a 
slow convergence of the expansion: c\/cq ~ —6.2, so the NLO term is as large as the LO for 
realistic values of a s . It has been argued [25,24] on the basis of various model calculations 
that the perturbative expansion of x can be improved, using nonperturbative information, 
in such a way that the perturbative expansion of A makes more sense. Here we prefer to 
treat A as an unknown free parameter. 

Since A is defined only through the formal all-order resummation eq. (3.7), its scale 
dependence is presumably non-perturbative. If xo is treated on the same footing as %o 
then, given that %o is scale-independent, A is effectively treated as being proportional to 
a s . In other words, even though the expansion AA eq. (3.4) starts at 0{a 2 ), AA is treated 
as being effectively order a s , so that xo is also scale independent. This approximation 
to the unknown non-perturbative scaling of A need not be correct. An alternative simple 
option consists of assuming that the value of A is scale-independent. Such an assumption 
can however only be valid as an approximation in a limited kinematical region, because 
asymptotic freedom implies that if we take the limit Q 2 — > oo at fixed x, then at some 
sufficiently large scale the low-order perturbative behaviour must be recovered. If A were 
strictly constant this requirement could only be satisfied if it also vanished: the 0(a s ) 
contribution to 7+ is ^yv^^x) ' anc ^ = ^ — a s^o which only vanishes as Q 2 — > 00 if A 
also vanishes. However, A approximately constant might be a a reasonable approximation 
over a limited range of Q 2 , and indeed in some explicit models [20], where the exact 
asymptotic small x behaviour can be computed, the value of A turns out to be reasonably 
well approximated by a constant [33] in the HERA region. We will thus consider both 
these options, and compare their phenomeno logical viability in sec. 5. 

The subtracted x% can then used to compute the resummed leading, next-to-leading,... 
singularities 7+, 7+, . . ., of the anomalous dimension 7+ by means of resummed versions 
of the usual duality relations [23,22]: instead of (2. 14), (2. 15), we now have 

Therefore, this resummation of the anomalous dimension amounts to a reorganization of 
the small x expansion of the anomalous dimension, i.e. the expansion of 7 (AT", a s ) in powers 
of a s at fixed : formally 

7 + (A, a.) = 7+d) + «-7+ (ft) + 0(a 2 s ) = 7+(ft) + a s ^+(^) + 0(a 2 ). (3.9) 

It is easy to see from (3.8) and (3.3) that to NLLx 

^+ (as \ _ ^gg ( qls \ 
is \ N ) Is y TV — AA J ' 



7ss ( TV ) 
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at LLx one simply lets N ^ N — AA in the unresummed anomalous dimension, while at 
NLLx one also lets xi ~ * Xi m (2.15). One can see explicitly that, because AA is formally 
of order a 2 ., the resummed and unresummed expansions eq. (3.9) are equivalent at NLLx, 
and differ only by terms of order NNLLx. 

3.2. Resummed structure functions 

Having discussed the resummation of the eigenvectors of the anomalous dimension ma- 
trix, we now turn to the determination of resummed structure functions. To this purpose, 
we need first to construct resummed expression for the individual parton distributions, and 
then resummed coefficient functions. In general, parton distributions can be determined by 
decomposing Q and G in terms of large and small eigenvector components, Q = Q + + Q~ , 
G = G+ + G~, with 

Q ± (N,t) = K^ g G ± (N,t), (3.11) 

which can be effectively done by means of the projectors M ± eq. (2.22). Because the 
coefficients K qg in general depend on t = InQ 2 , they also contribute to the scale dependence 
of the parton distributions. The coefficients K^ g were given to leading nontrivial order in 
eqs. (2.4), and are explicitly determined using the color-charge relation eq. (2.2) and the 
quark sector anomalous dimensions, given in the DIS scheme in eq. (2.18). In fact, to 
NLLx eq. (2.4) implies that 

j f G+(x, t) = 7 ^ LLx (iV, a s (t))G+(x, t) + 0(a 2 s ), (3.12) 

so the determination of the quark and gluon distribution is straightforward. 

Now, the quark anomalous dimensions ~f qg , ^ qq eq. (2.18), and the coefficient K+ g 
eq. (3.11), and indeed the longitudinal coefficient functions eq. (2.27), as well as all the MS 
coefficient functions, are all determined [29] as functionals of the leading log x anomalous 
dimension 7+ (jf)- It is easy to understand the reason for this as a further consequence of 
the kr factorization [29] which gave the duality [19,23,22] equations (2. 14), (2. 15) relating 
7 + to x- We thus quickly review this derivation [23], after which the natural extension of 
our resummation of 7 + to the quark sector will become clear. At small x and large Q 2 , kr 
factorization implies that we can write the large component Q + of the quark distribution 
in the factorized form 

Q + (N,t)= / ^—e Mt K+(M,N)G+(M,N), (3.13) 
where the Mellin transform has been defined by 

/oo 

dte- Mt Q + (N,t), (3.14) 
-00 

with inverse 

Q + (N,t)= ^e Mt Q+(M,N), (3.15) 

Jc—ioo 
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where the contour passes to the right of any perturbative singularities near M = 0, but 
to the left of M = |. The leading twist behaviour of parton distributions is then given by 
closing the contour to the left and picking up the contribution of the rightmost singularity 
of the integrand in eq. (3.15) (see e.g. refs. [19,22]). The large component of the gluon 
distribution G + (M, N) satisfies the BFKL equation, with kernel %(M, a s ) eq. (3.1) and 
LLx solution 

G+(M, N) = - Gq+(M Lv (3-16) 
V ' N-a sXo (M) V ' 

It thus has a simple pole at M = 7+ (AT, a s ), where 7 + is given by the duality relation 
eq. (2.14). Solving the BFKL equation at higher orders gives corrections to the location of 
the pole as a series in a s at fixed which correspond to the duality relations eq. (2.15), 
whereby all anomalous dimensions are expressed as functions of the position of the LLx 
pole. 

Because K+ g (M, N) eq. (3.13) is free of collinear singularities [29], and is therefore 
regular in the neigbourhood of M = 0, the same pole which gives the leading twist small 
x behaviour of G + (N,t) also gives that of Q + (N,t). It follows that 



Q+(N, t) = K+(>y + (N, a s ), N)G+(N, t) 

= ^(7 s + (f))G + (jV,f) + 0(a s ), 

where we have replaced 7+(iV,a s ) by 7+(ff), and K+(M,N) by K+(M,0) = K+(M) 
since the explicit N dependence only generates subleading singularities. Thus the coeffi- 
cient K+ g at leading nontrivial order only depends on through 7+. Explicit computation 
of the cross-section for off-shell photon-gluon scattering [29] in MS and DIS schemes gives 
(in DIS) K+ g (M) = h 2 (M)R(M)/M, which when combined with (3.17) lead to the ex- 
pression (2.18) for the quark anomalous dimension 7|| . A similar argument shows that all 
the leading MS coefficient functions and anomalous dimensions have the same property, 
i.e. they may be determined as functions of M which is then at leading twist set equal to 
1s~(~n) by the duality relation at LLx. 

After the resummation eq. (3.9), the LLx anomalous dimension becomes 7+ eq. (3.10). 
The resummed expression of parton distributions to NLLx are then simply found by using 
the resummed anomalous dimension (3.10) not only for the evolution of the large eigen- 
vectors, but also in the expression eq. (2.18) for the quark anomalous dimension, and thus 
the coefficient K+ g eq. (3.11), (2.4). The resummed quark anomalous dimensions are thus 
given by 

«,7 s 9 f m = ^(7 s + (t))^(7 s + (^)) = a s ^J (^), (3.18) 

since the pole in (3.16) shifts from 7^(7^-) to 7^(7^) = 7^~( tv-aa ) - Thus the resummation 
in the quark sector leads simply to a shift N ^ N — AA in the argument of the quark 
anomalous dimensions. 

The resummed expression of the structure function F2 in the DIS scheme is given by 
the resummed quark distribution. Likewise, we can obtain a resummed expression for the 
longitudinal structure function, by using kr factorization to express, at small x and large 
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Q 2 to NLLx the large component of F L (N,t) as 



rc+ioo 
J c— too 



dM 
2^i 



,Mt 



(KlJM)Q+(M,N) + Kl as (M)G+(M,N)), (3.19) 



where the coefficient functions are free of collinear poles and given by [29] SS (M) = 
hL(M)R(M) and a color-charge relation. Taking the residue of the pole then gives (2.27) 
for the coefficient functions in the unresummed case, whereas after resummation 



N—A\ 



)• 



(3.20) 



Thus the resummation again amounts to a shift N — > N — AA in the argument of the 
coefficient function. Clearly the same holds true for MS F2 coefficient functions, and 
indeed for any hard cross-section which is free of collinear singularities. 

3.3. Double leading resummation 

Having constructed resummed expressions for anomalous dimensions, coefficient func- 
tions and structure functions to NLLx, we can now combine them with the two loop results 
according to the lines discussed in the previous section to construct a resummed double- 
leading expansion. As discussed in ref. [22] , the resummed double-leading expansion of the 
large eigenvalue 7^ L R is constructed by replacing the singular terms 7 S , 7 SS ,. . . in eq. (2.6) 
with their resummed expressions eq. (3.9), i.e., in practice, by letting N — > N — A A in the 
corresponding contributions: 



a sl +(N) +1 + (jfin 



+ a s 



- a. 



N 



- e r 



(3.21) 



+ 



Likewise, we can construct the quark sector double-leading anomalous dimensions and 
coefficient functions by performing the same replacement in eqs. (2.9), (2.26): for example 
the resummed double leading anomalous dimension 7^ r is given at NLO by 



and similarly for 7^ R , while 



_ e£ q 

ss \ N-AX I as ~N e 



+ 



CI, DL)R (iV, a s ) = a s [CI X {N) + ss (^) - e\ 



(3.22) 



(3.23) 



Notice that the momentum subtractions eq. (2.12), (2.13) are affected by the resummation 
and therefore must be recomputed at the resummed level. 

As already discussed in ref. [22], this procedure generates an ambiguity in the treat- 
ment of the double-counting terms in 7d L : because these terms are common to the loop 
expansion 70, 71,. . .and the small x expansion 7 S , 7 SS ,. . . , we are free to decide whether 
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to leave them unaffected by the replacement N ^ N — AX (as in eq. (2.6)) or not. The 
difference between the two procedures is formally sub-subleading, provided we suitably 
readjust the subleading double-counting subtractions. An equally acceptable alternative 
to the 'R-resummation' eq. (3.21) is thus the 'S-resummation' 



Oi 8 ll(N) + 1ss( N 1 s AX ) + Tr(jV-AA) 2 ^ s N — AX e Q 



(3.24) 

+ •••, 



where the double-counting subtraction is now also resummed. Clearly, a variety of inter- 
mediate alternatives would also be possible. If the S-prescription eq. (3.24) is used, then 
the double counting terms in the quark sector anomalous dimensions are also affected by 
the replacement N — > N — AA. Since in these terms all singularities start at the NLLx 
level, no further readjustment is necessary: 

7£i,s W «.) = ^ q 9 (N) + a s [a sl r(N) + 7 ™ (^) - - eg] + • • • (3.25) 

The resummed longitudinal coefficient functions (3.23) are unchanged in the S- 
resummation, since there the double counting terms are independent of N. 

The main difference between the two resummed expressions eq. (3.21) and (3.24) is the 
nature of the small-a; singularities of 7+: the resummed anomalous dimension always has 
a cut starting at N = A, which corresponds [23] to an x~ x behaviour of splitting functions 
at small x. However, if the S-prescription eq. (3.24) is adopted, the anomalous dimension 
also has a simple pole at iV = 0, which leads to a "double-scaling" [14,15] rise at small x. 
If A is positive, then the power rise will dominate the asymptotic behaviour, and the two 
resummations give similar results [22], but if A < the double scaling rise is dominant. In 
the latter case, the S-resummation eq. (3.21) will give results at small x which are very 
close to those obtained in a fixed leading or next-to-leading order computation, which are 
dominated by double scaling, unlike the R-resummation eq. (3.24), which instead would 
give a valencelike drop of the splitting functions at small x. Indeed, even when A = the 
R-resummed splitting function and structure function drop logarithmically [23]. 

It turns out that different prescriptions for the all-order running of the coupling, to 
be discussed in sec. 4 below, can lead to anomalous dimensions whose small N behaviour 
is characterized by either a pole or a cut; we can therefore take our R vs. S ambiguity 
in the resummed NLLx double-leading result as an indication of this further ambiguity, 
which could only be resolved by arguments that go beyond perturbation theory. The 
phenomenological implications of these different options will all be discussed in detail in 
sec. 5. 



4. Scheme dependence and resummation ambiguities 

In this section, we discuss the dependence of the resummed structure function on the 
choice of factorization scheme and its interplay with the resummation procedure. At the 
double-leading level, a factorization prescription must be specified both in the treatment 
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of the leading InQ 2 and of the leading Inx terms. The first choice allows one to switch, 
for instance, from the DIS to the MS scheme, while the latter choice relates different 
small x resummations such as DIS and QoDIS [34]. This latter choice interferes with 
the resummation prescription, thereby raising theoretical issues such as the uniqueness of 
the resummation procedure, and the process-independence of the resummed results. This 
will also lead us to discuss some ambiguities which are related to the way the running 
of the coupling is treated in the Regge limit. These issues are not directly relevant for 
phenomenology, but they are useful in order to clarify the theoretical underpinnings and 
limitations of our approach. Readers who are only interested in the construction of the 
standard MS scheme can skip this somewhat technical discussion (contained in sec. 4.2- 
4.4), while those who are interested only in the resummed DIS results should skip this 
section altogether. 

4-1. NLO scheme changes 

In general, to next-to-leading order in the double-leading expansion, we can define a 
double-leading scheme change matrix eq. (2.25) as (in matrix notation) 



U(N,a 8 ) = 1 +a a z DL (N,a a ), z DL (N, a s ) = Zl (N) + z ss (%) - z°, (4.1) 
where Zl (N) is regular at iV = 0, z ss (jf) = E7=o z ss iwY '» and z i = Z M = z ss The 



change in the double-leading anomalous dimension matrix (2.8) induced by using this form 
of the scheme change matrix in eq. (2.25) is then 



with a corresponding change in the subleading double counting term. The scheme change 
separates into a standard NLO loop scheme change, and a small x scheme change [28] 
(z ss being the V matrix of ref. [28]) because all cross terms are formally subleading: 
while the constant contribution z® contributes both to the 71 and 7 SS scheme changes, 
z\{N) — Zi introduces extra factors of N and thus cannot contribute to a change in 7 SS 
while z ss (jf) — zi introduces extra factors of a s and thus cannot contribute to a change in 
71. Furthermore, the (3q term only contributes to the two-loop part of the scheme change 
eq. (4.2), because the running of the coupling is a NLLx effect. 

Let us now consider the specific case of the DIS— >MS scheme change. The MS scheme 
is defined by the computation of the corresponding collinear finite coefficient functions, 
whose singular terms have been determined in ref. [29] . The F 2 double-leading coefficient 
functions are then 



where i = q, g, and there is a color-charge relation C^, ss = §j(Cl ss ~ ef) just as in the 
longitudinal case, eq. (2.26). The double-counting subtractions are now ef = e\ = 0, 



I'l = 7i + [*i,7o] - 

l'ss = Iss + [Zss,ls], 




(4.2) 



CI DL (iV, a.) = 5* + a s [C\ X (N) + C\ ss (f ) - e\] , 



(4.3) 
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and the coefficient functions may be evaluated through a series expansion analogous to 
(2.28): 

oo 

^ 2 %(7f) = ^(l + £^ S (fr), (4.4) 

n=l 

where the coefficients d^ s are again listed in table 2. Unlike the longitudinal coefficient 
functions, which start at next-to-leading order in the double-leading expansion, and there- 
fore are unaffected by scheme change at this order, the F-i coefficient functions begin at 
leading order and thus depend on the choice of scheme. 

Demanding invariance of F<i implies that the scheme change matrix which takes us 
from the DIS scheme eq. (2.24) in which DL = $ iq to the MS scheme in which eq. (4.3) 
holds must have 

U q9 = -Ci» L tf™ = 2-C* DL . (4.5) 

Because the Fi coefficient functions have the form eq. (4.3), and in particular they do 
not contain any leading singularity, i.e. C2, s = 0, the quark entries of the scheme change 
have the form of eq. (4.1). These entries are sufficient to determine [28] the singular 
contributions to the MS quark anomalous dimensions eq. (2.8): 

i q s 9 s ' = il! + « 5 , ilf = m + fr--f- ( 4 - 6 ) 

These are the same as the leading singularities of the MS quark anomalous dimension 
computed in ref. [29]. In other words within the double leading prescription the quark 
entries of the double-leading anomalous dimension matrix eq. (2.9) transform consistently: 
the DIS and MS double-leading expressions are found by combining the respective DIS 
and MS large x and singular terms. 

Since the MS quark anomalous dimensions still satisfy the color-charge relation (2.5), 
it follows that the singular terms in the small eigenvalue remain zero in MS. Moreover 

„99' _ ~gg _ Cp_ gg gg (a y\ 

Iss Iss Ca ss ' s ' ^ ' 

so it follows from eq. (2.3) that the singular part of the large eigenvalue remains unchanged 
too. This latter conclusion may also be seen directly from (4.2): since the scheme change of 
7 S + a s 7ss is a commutator, whose trace vanishes, the sum of the eigenvalues must remain 
unchanged, so 7+ + 7" = 7+ + 7^. Since the singular parts of 7" remain unchanged, 
those of 7 + must too. 

In order for the full scheme change to be consistent, it is sufficient to require that 7^ 
coincides with the usual MS matrix of NLO anomalous dimensions, and that the scheme 
change matrix be of the form eq. (4.1). Indeed, it then follows that the eigenvalues of 71 
transform in the standard way, while the eigenvalues of 7 SS do not change at all. Because, as 
discussed in sec. 2, the eigenvalues of the double-leading Y DL are the sum of the eigenvalues 
of their small x and large x components, this ensures that the double-leading eigenvalues 
eq. (2. 6), (2. 7) also transform consistently, i.e. they are the sum of the large and small x 
contributions in the corresponding scheme. This is sufficient to fix the scheme, since, given 
eigenvalues and quark entries the evolution equations and their solutions are completely 
determined. 
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The remaining two entries of the scheme change matrix may be fixed entirely by a 
natural extension of the method employed at two loops [32] : this gives 



U 9q = Cl^-l tf" = l + Cf >DL . (4.* 
In fact the only effect of these entries beyond two loops is to change 7f|: 



r^m' — rydl I ry99 
ISS ISS • IS 



£f_ ( 7 99 _ 7 99\ _ 7 99 (A Q) 

C A ^ ss ss ' ss ' \^-^J 



However, as we explained in sec. 2, 7f| is of no consequence at NLLx. 

4-2. LLx scheme changes 

On top of the 'standard' double leading scheme changes eq. (4.1), if we allow any 
factorization scheme change of the form eq. (2.25) provided only that the leading log x 
anomalous dimensions be unchanged, there is the freedom to perform an extra 'leading 
log x } scheme change. This is a consequence of the fact that the leading small x singularities 
only appear in the gluon sector of the anomalous dimension matrix eq. (2.1). Indeed, this 
requirement is satisfied [28] by any scheme change of the form 

U= (^/(a 8 /N) zf(aV))' (4 ' 10) 

where z 99 (a s /N) = 1 + z 99 1 jf + . . . . This scheme change amounts to a LLx modification 
of the normalization of the gluon distribution, so the the DIS identification of F 2 with the 
quark distribution remains unaffected. 

Upon this scheme change 7" and 7" 1 " change according to 



I q s 9 s' = i q s 9 s/U{lt) 



Tss Tss ~r 



^X'oils) dM 



(4.11) 



where we have for convenience defined u(M) by the implicit equation z 9 s 9 '(a s /N) = 
u (ls(jf))j an d made use of the fact that (differentiating (2.14)) 



51n « s x'odl) 

The anomalous dimension 7|| changes in such a way as to preserve the colour-charge 
relation (2.5), and thus 7" is unaltered. The scheme change of the NLLx anomalous 
dimension eq. (4.11) corresponds to a dual scheme change [23] of the next-to-leading BFKL 
function 

x[(M) = Xi(M)-^ X o(M) d ^. (4.13) 

Since these scheme changes only affect the singular contributions to the anomalous 
dimensions, for the remainder of this section we will only discuss the small x expansion 
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of the structure function. In fact, they generate a class of schemes which only differ in 
the way the small x resummation is f actor ized, and specifically, in the factorization of 
the expression eq. (3.13) of Q + in terms of the solution to the large eigenvector evolution 
equation and the coefficient K+ g . For definiteness, let us assume that we start from the 
DIS scheme, where F 2 and Q coincide: the LLx scheme change then moves within the class 
of DIS schemes. 

Upon resummation in any of these schemes, the expansion of x gets rearranged accord- 
ing to eq. (3.2), and the subleading anomalous dimensions eq. (3.10) change accordingly. 
However, the scheme changes do not in general commute with the resummation procedure, 
in the sense that a scheme change characterized by a given function u(M) does not have 
the same effect if performed before or after resummation. To see this, consider a pair of 
schemes connected by a scheme change function u(M), whose NLLx anomalous dimensions 
are thus related by eq. (4.11). We can now determine the resummed anomalous dimen- 
sion in the primed scheme in two different ways. One possibility is to first perform the 
resummation eq. (3.10) in the unprimed scheme, then transform the result according to 
eq. (4.11), thereby obtaining a NLLx anomalous dimension [7+]': 



[its]' = 



Xi (7+) , A) 1 d[a sXo (M)} dlnu 



+ 



Xo(7* + ) 47r Xo(7s + ) da * 



dM 



M- 



=7+' 



(4.14) 



where Xo(M) eq. (3.3) depends on a s through AA. 

Alternatively, we first perform the scheme change eq. (4.11) of the unresummed 
anomalous dimension, and then resum the result according to eq. (3.10), thereby obtaining 

a NLLx anomalous dimension (stable in the sense of eq. (3.5)) [ 7a " 



[7, 



ss ] 



Xi(7 



1 



Xo(7.) 4 ^Xo(7 S + ) 



Xo(7 s + ) 



dlnu 



dM 



M=7, 



_ + -xo(|) 



dlnu 



dM 



(4.15) 



The results obtained by the two procedures (4.14) and (4.15) thus differ by 



[its'] - [its]' = 



A) 1 
47r Xo(7 



^ Xo(7* + ) 



dAX dlnu 



da q dM 



M-- 



dAX ( dlnu 



da* \ dM 



=7 
dlnu 



. 1 din u 
+ Xo( 3 ) W - 



m=y s 



dM 



M=\ 



+ 



d\ dlnu 



da? dM 



M=: 



(4.16) 



Because A is formally of order a s and A A of order a^, the non-commutativity given by the 
first term on the r.h.s. of eq. (4.16) is formally subleading, whereas in the second term it 
is leading (i.e. the same order as 7^). 

To understand this, notice that at 7^" = | the second term in eq. (4.14) will in general 
be singular, so [7^]' eq. (4.14) will be unstable, unless 



dX dlnu 



da, dM 



= 0. 



M=- 



(4.17) 
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By contrast, [7^'] eq. (4.15) is stable by construction. This means that if we have re- 
summed the anomalous dimension to give a stable expansion, after a further scheme 
change the expansion will still be stable only if the condition eq. (4.17) is satisfied. Oth- 
erwise stated, the subtraction coefficient c\ eq. (3.4) required to stabilize the primed and 
unprimed schemes is not the same. The difference between these stable and unstable 
prescriptions is the leading order second term on the r.h.s. of eq. (4.16). The condition 
(4.17) for the vanishing of this term may thus be seen as a condition on the further scheme 
changes which are permissible after resummation has been performed. 

After resummation, as discussed in sec. 3.1, the scale dependence of A is unknown: 
we will be model it by taking A to be either proportional to a 8 , or scale-independent, even 
though, as discussed above, the latter option can only be true as an approximation valid in a 
limited kinematical range, or if A = identically. In both cases AA contains a contribution 
of order a s : if A oc a s then ^ = while if 4^- = 0, then = _ A<i . Note that this 

* * da s ol s ' aoc s ' doL a a s 

implies that the size of this contribution is comparable in the two cases whenever |A| <C Ao: 
this will turn out to be the case in actual phenomenology (sec. 5). In both cases, the first 
term on the r.h.s. of eq. (4.16) leads thus to a potentially significant ambiguity in the 
scale dependence of resummed parton distributions upon LLx scheme changes eq. (4.10). 
We have however checked that for common scheme choices this ambiguity is significantly 
smaller than that which we explore by switching from the S- to the R-resummations 
discussed in sect. 3.2, and concentrated in the same region. On the other hand, the second 
term on the r.h.s. of eq. (4.16) is proportional to -fj^- = if A oc a s , while it vanishes 
if A is independent of a s (and thus in particular if A = 0). Hence, the stability condition 
eq. (4.17) is automatically satisfied in any scheme if A doesn't depend on a s . Furthermore, 
the condition is approximately satisfied even when A is proportional to a s , provided it is 
small enough: in such case, the ratio eq. (3.5) rises linearly P S8 (x, a s )/P s (x, a s ) ~ In 1/x, 

but the slope of the rise is proportional to A [23] . 

4-3. Coefficient function resummation 

So far we discussed the implications of LLx scheme changes for the stability of the 
resummed process-independent large anomalous dimension 7+ in DIS (or MS). However, 
as discussed in Sect. 3.2, the coefficient K gg eq. (3.11) and the coefficient functions also 
contribute to the scale-dependence of the structure function. In order to achieve pertur- 
bative stability of the structure function, therefore, we must consider the resummation of 
these process-dependent quantities as well. To do this, we first define an 'effective' NLLx 
anomalous dimension 

7iSLL*,affW «-(*)) = f t ^F+(x,Q 2 ), (4.18) 

where for definiteness we consider first the structure function = Q + in DIS schemes. 
The effective anomalous dimension is by definition scheme-independent, but process de- 
pendent: it coincides with 7^ LL x on ly m schemes where not only the coefficient functions 
are trivial (such as DIS), but also the coefficient K+ g eq. (2.18) is trivial, in the sense that 
neither contributes to the scale dependence of F 2 (x, Q 2 ). It follows from eqs. (3.12), (3.11) 
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that we can always go from the resummed DIS anomalous dimension to the effective 
anomalous dimension by means of a LLx scheme change (4.11), with 



' a,M X0 (M) ' 



where we have used the fact that a s Xo(7^) = N independent of t, and the normalization 
constant n qg = (with e\ defined in eq. (2.9)) is to ensure that the scheme change takes 

the required form (4.10), which requires u(0) = 1 [28]. After resummation in DIS 



TNLLx.effC^"-) =7NLL x ( iV 5 " S ) + 

1 (4-20) 



for 1+ 1 d(a s x )d\n(h 2 (M)R(M)/M) 
4tt I 



^(7+) da s dM 



At 7+ = | the second term is singular, and the effective anomalous dimension unstable, 
unless 



dX 

da« 



n d\nR-\ 
2 H 



= (4.21) 



in the limit M — > |. Notice that, because h 2 {M)/M 2 is symmetrical about M = |, the 
stability condition (4.21) only involves the process independent function R(M). Similar 
results may be derived for the unstable contribution of the coefficient function to the scale 
dependence of other physical observables: that for F L turns out to be the same as (4.21), 
whereas for other quantities, such as for instance heavy quark production cross-sections, 
they would be different. 

One might think therefore that the effective anomalous dimension (and thus the re- 
summed expression for, say F 2 ) may be stabilized by resumming the process-dependent 
coefficients as well, i.e. by a further subtraction analogous to c\ eq. (3.3). If we try to do 
that, however, we have a problem: R(M) ~ (| — M) -1 / 2 when M ~ \ [29]? so the required 
subtraction diverges in the M — > \ limit. This is to be contrasted to the case discussed 
in sect. 4.2, eq. (4.16), where even though [7^]' was singular in the limit M — > ^ the 
subtraction needed to stabilize it is finite. 

There are several ways of dealing with this issue: 

(i) We may argue that it should be ignored, on the grounds that the process-dependent 
F2 coefficient is actually leading: after all, without it F 2 would vanish. Hence, the effective 
anomalous dimension argument is not relevant: we have two contributions to scaling vio- 
lations of physical observables, one from the evolution of the parton distributions (which 
should be stable) , but another from the hard cross-section, which may be relatively large. 
Stability of the latter becomes a relevant issue only when the next-to-leading correction to 
it is computed [35], included in it (at the next order), and compared to the leading order. 

(ii) We may ignore it, on the grounds that the effect is in practice rather small. 
Indeed, if A is taken to be proportional to a s , when the stability criterion eq. (4.21) is 
not met, so for the effective splitting function P s e f (x, a s )/P s eff (x, a s ) ~ (lnl/x) 3 / 2 [23] 

x— >0 

the slope of the rise is proportional to A. Now, it turns out (sec. 5) that the actual 
phenomeno logical value of A is generally close to zero, so the slope of the rise is rather 
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small. Of course, the rise goes away altogether if A is taken to be independent of a s . In 
fact, one can check explicitly that this rise is undetectable for all practical purposes in the 
HERA region provided |A| ^ 0.2. Therefore, for phenomenological applications we may 
stick to the resummation procedure discussed in sec. 3, and simply check a posteriori that 
no significant instabilities are present. 

(iii) If A = 0, it makes no difference whether the subtraction is determined from 
the effective anomalous dimension or in MS, since then the stability criterion (4.21) is 
always met. Various theoretical arguments [19] suggest that A = 0. Furthermore, in the 
S-resummation, A ~ is consistent with the HERA data, as we will see below. 

(iv) One should trace the origin and meaning of the singularity in -R(M), and try to 
deal with it through a further resummation. This will be discussed in sect. 4.4, where we 
will show that the source of the singularity appears to be related to the NLLx truncation of 
running coupling effects in the BFKL kernel. Therefore, a model-independent resummation 
of the singularity does not seem possible at present, though we will argue that some effects 
of it can be modeled by comparing the S- and R-resummations of sec. 3.3. 

4-4- Regge limit and the running of the coupling 

The singularity in the process-independent function R(M) is related to some limi- 
tations of the NLLx approach, and thus warrants further discussion. The origin of this 
singularity is exposed by noting that R(M) = r(M)/ a/— x'o(M) where the function r(M) 
is regular if < M < 1. Then, the function Xeff which is dual to the effective anomalous 
dimension (4.18) contains a term (see eq. (4.13)) 

This contribution to Xeft has a simple pole at M = |. As discussed in sec. 4.3, the 
presence of a pole in xi at M = \ is a problem in view of the fact that the Xi{\) fi xes the 
value of the NLO subtraction c\ eq. (3.6) which would be required to stabilize XeS- This 
problem will be present in schemes (such as the QoDIS scheme [34]) where the contribution 
eq. (4.22) is included in the anomalous dimension, and indeed in any scheme if we attempt 
to implement a resummation prescription at the level of effective anomalous dimensions, 
as discussed above. 

The origin of such a singular contribution to x can be traced to the perturbative 
treatment of the running of the coupling to NLLx. Indeed, the effects of the running 
of the coupling can be viewed as the addition of higher order terms to the anomalous 
dimension. It turns out [23,22] that these terms can be always expressed as an order- 
by-order modification of the duality relations: the form of the duality relation remains 
unchanged, but the function %(M) no longer coincides with the BFKL kernel, rather, it 
differs by it by the addition of contributions of order ((3oa s ) k . To NLLx, there is thus 
a running coupling 0(/?o) contribution to xii which turns out [23] to be precisely of the 
form (4.22). The singularity can be moved from the anomalous dimension to the coefficient 
by a LLx scheme change but it is of course always present in the effective anomalous 
dimension. In fact, this is what happens in the DIS (or MS) scheme, where the singularity 
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in the running coupling contribution to Xi is cancelled by a similar contribution from 
the normalization of the LLx unintegrated gluon distribution, but then re-appears as the 
square root singularity in the R(M) contribution to the coefficient function K qg eq. (3.17). 
All this suggests that the running coupling singularity eq. (4.22) and the singularity in the 
DIS stability condition eq. (4.21) are one and the same. 

The running coupling singularity eq. (4.22) signals a failure of the NLLx treatment 
of the running of the coupling in the vicinity of the point M = |. In this respect, it 
is analogous to the singularity in the NLLx anomalous dimension 7^(7^) eq. (2.15) as 
7+^| which leads [23] to the instability which is removed by the resummation of x 
eq. (3.2). This analogy suggests that the running coupling singularity (4.22) might be 
understood, and possibly removed, by treating the running of the coupling to all orders. 
Although interesting models for the running of the coupling in the BFKL equation beyond 
NLLx have been discussed [24], a model-independent treatment does not appear to be 
possible at present. An investigation of models is beyond the scope of this paper, where 
we only wish to make use, as far as possible, of model-independent results. However, it is 
worth noticing that the inclusion of running coupling terms beyond NLLx leads to a series 
of contributions to the duality relation to all orders in /3o, of which x rc eq. (4.22) is the 
leading term. Resumming this series to all orders may change the qualitative structure of 
X in the vicinity of the leading order minimum. 

For instance, if we assume that we can let the coupling run by setting 

a s -> ^ — (4.23) 

1 "I" 47r «s dM 

it is possible to determine in closed form [36] the solution G(M, N) to the BFKL equation. 
If one further approximates the kernel x with a parabolic form x = X + k(M — |) 2 , which is 
accurate close enough to the minimum, it is possible to determine analytically the solution 
G(N,t) in terms of the Airy function [37], and thus the associated anomalous dimension. 
One can then show that whenever A > 0, the anomalous dimension has a simple pole in 
the range < A < A. This changes qualitatively the perturbative behaviour of the small x 
terms 7 S , 7 SS ,. . . : for all A > 0, if x is parabolic then 7 which is dual to it has a branch point 
at N — A, however the branch point turns into a simple pole when running coupling effect 
are resummed according to eq. (4.23). Now, recall that the simple one-loop anomalous 
dimension is characterized by a simple pole at N = 0. I follows that if A is close enough 
to zero, the small x behaviour induced by this resummation of running coupling effects, 
which also generates a pole in the anomalous dimension, will be quite close to that of the 
one-loop result. This is also the small x behaviour which is found in the S-double leading 
expansion eq. (3.24), while in the R-double leading expansion, as discussed at the end of 
sec. 3, there is no pole but only a cut. Therefore, switching between the two resummed 
forms of the double-leading expansion discussed in sec. 3 changes the form of the dominant 
small N singularity in a way which can give a feeling for the kind of uncertainty related 
to the resummation of running coupling effects in the Regge limit. 

5. Comparison of the resummed structure functions with the data 

In the previous sections we have developed a formalism to resum a sequence of all 
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order small x corrections to splitting functions and coefficient functions. The resulting 
expressions for structure functions should be valid in an extended range of x down to 
much smaller values than would be the case in the usual two loop approximation. We now 
compare our theoretical description of structure functions with the data, with particular 
attention to the HERA data which probe the small x region. When compared to the two 
loop results, the resummed calculation has several important new features. First, in our 
formalism, in addition to the input parton densities at some initial scale Qq and the value 
of a s (m 2 ), we also have to extract from the data the value of the exponent A that fixes the 
asymptotic behaviour of the resummed splitting functions in the Regge limit. Second, our 
resummation procedure has a number of further ambiguities which must either be resolved 
on phenomeno logical grounds, or else considered to be part of the theoretical errors on the 
fitted quantities. Specifically, as discussed in sec. 3, there is some uncertainty in deciding 
the scale dependence of A. Furthermore, there is an important ambiguity in the treatment 
of the double-counting terms in the resummed double-leading expansion, which gives rise 
to the pair of alternative prescriptions, denoted by R and S, and discussed at the end of 
sec. 3, eqs. (3.21), (3.24). 

Since this work is focussed on small x data, we want to start from input parton 
densities that already embody the information contained in the large x data so that we can 
study the details of the small x behaviour without having to perform a global fit. Therefore, 
our procedure for analysing the small x data is to start from a set of globally fitted parton 
distributions, at a starting scale Qq, and treat as free parameters the exponents \ q and X g 
which characterize the small x behaviour of the singlet quark and gluon as Q = xq ~ x~ Xq , 
G = xg ~ x~ Xg . Of course for each new choice of X q and X g all the other parameters 
must be readjusted in order to maintain the large x shape (say for x ~ 0.01) of the parton 
distributions and preserve the momentum sum rule. These new distributions, together 
with the valence and nonsinglet densities, which are left unchanged, are then used as input 
distributions, evolved up to the HERA data, and X q and X g are tuned to obtain a best 
fit. Of course, X q and A s should not be confused with A: the former are parameters in the 
input parton distributions at Qq, while the latter is a property of the resummed splitting 
functions and thus of the parton evolution. 

The data set that we use consists of HI data collected in 1995-1997 [38]. The HI data 
have been recently shown [39] to be consistent with large-ir data as well as with the LEP 
value of a s , which dominates the world average [40]. We only include neutral current data 
with x < 0.01 and 5 GeV 2 < Q 2 < 1000 GeV 2 in order to be above the charm threshold 
and well below the Z. We fit directly to the 'reduced cross-section' 

a red (x,y,Q 2 ) = F 2 (x,Q 2 ) + V \ 2 F L (x,Q 2 ). (5.1) 

2(1 -y) + y z 

The x 2 is calculated with statistical and systematic errors added in quadrature, except for 
the overall correlated normalization uncertainties which are treated separately. As input 
partons we use MRSA4 [41], which have ot s {m 2 ) = 0.120. We checked that essentially the 
same results are obtained with the CTEQ4A4 set [42] . 

5.1. Fixed a s {m 2 ) 

In a first round of fits we fix a s {m 2 z ) = 0.119, a value close to the current world 
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100 




Figure 1: \ 2 f° r the fit to the 95 HI data [38] of the reduced cross sec- 
tion eq. (6.1) at two loops (solid), and in the double-leading expansion, S- 
resummation (dotdash) and R-resummation (dashes) as a function of A. Here 
a s {m 2 z ) = 0.119. 



average [40], and proceed to fit \ q , g at Qo = 2 GeV for different values of A. The number 
of degrees of freedom is in this case ndof = 93. In fig. 1 we present the behaviour of 
the best-fit % 2 as a function of A for different variants of our approach. In fact, we show 
curves for x 2 computed in the two resummed expansions discussed at the end of sec. 3, 
as a function of A, and compared with the two loop value. In the resummed fits, we have 
taken X(Q 2 ) = a s (Q 2 )c, where c is a constant independent of Q 2 . The value given on the 
abscissa in fig. 1 is (for definiteness) A = 0.2c; 0.2 is taken as a representative value of 
a s (Q 2 ) in the HERA data region. 

We see that in the case of the S-resummation a fit at least as good as the NLO fit 
is obtained over a wide range of negative values of A. In particular the resummed fit has 
the lowest \ 2 for A ~ —0.25. On the contrary the R-resummation only approaches the 
NLO result in a very narrow range around A w 0.2. Thus, it turns out that fitting the 
data with this resummation prescription requires a substantial fine tuning in A, which is 
however not necessary when using the S-resummation. This can be understood on the 
basis of the results of ref. [22], where it was shown that the S-resummation gives a small 
x splitting function which is extremely close to the two loop one whenever A < 0, whereas 
the R-resummation only gives a result close to the two loops one for a fine-tuned value 
of A 0.2. Comparison with the data indicates that only very small deviations from two 
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Figure 2: Starting gluon slope, G{x,4 GeV 2 ) ~ x~ Xg for the fit to the 95 HI 
data [38] of the reduced cross section eq. (6.1) at two loops (solid), and in 
the double-leading expansion, S-resummation (dotdash) and R-resummation 
(dashes) as a function of A. Here a s {m 2 z ) = 0.119. 

loops are tolerable. 

The behaviour of the best-fit gluon exponent X g at the initial scale Q = Qo is shown 
in fig. 2 as a function of A. In general, X g increases with decreasing A: the input gluon 
becomes more valencelike if the dynamical exponent from the evolution becomes larger. 
The fall of X g in the R-resummation is extremely steep: the input gluon changes rapidly 
in an attempt to compensate for the variation with A of the singlet evolution. This is again 
a sign that this expansion requires fine tuning. However, for the fine tuned best-fit A, the 
resummed value of X g is close to that obtained in the two loop fit. The results for X g found 
in the S-resummation are much more stable. In particular, we find that the input density 
is rather more valence-like than the unresummed two loop one throughout the region of 
small A where the x 2 is near the minimum. In fact, all the fits with a reasonably good % 2 
have a value of X g such that the input gluon density is flat or valencelike, and thus leads 
to double scaling behaviour (faster than any log but slower than any power) [14,15] when 
evolved into the HERA region. 

Note that at negative values of A the S-resummation inevitably approaches the two 
loop fit. This is because the resummation suppresses the high order terms when A — > — oo, 
so that the unmodified one and two loop terms dominate at small x, and the two loop 
results are recovered. Of course, we do not consider very large negative values of A to be 
realistic: the first and second terms of the series for A suggest that its value should be 
perhaps A <^ —0.5. Finally, note that at large values of A, the S-resummation and the 
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X 

Figure 3: The gluon distribution and longitudinal structure function corre- 
sponding to the best-fit value of X g of fig. 2 at Q 2 = 4, 20 and 100 GeV 2 at 
two loops (solid), and in the double-leading expansion, S-resummation (dot- 
dash) and R-resummation (dashes). 
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R-resummation come together because the more singular behaviour of the higher order 
resummed terms makes the ambiguous terms negligible at small x, as discussed in the end 
of sec. 3. 

The resummed gluon density and longitudinal structure function obtained in these fits 
are displayed in Fig. 3a,b as functions of x for three different values of Q 2 . In particular, we 
compare results obtained in the two loop fit with the best fits in the S- and R-resummation. 
We see that even though at the starting scale Qq = 4 GeV 2 both the resummed gluon 
densities are more valencelike than the two loop one, when evolved to higher Q 2 the R- 
resummation quickly becomes more singular. The same features can be seen reflected in 
Fl- This shows that the gluon extracted from a conventional unresummed fit could be 
sizeably different from the resummed one, which is sensitive to the small x corrections to 
the evolution equations. It also suggests that more precise data on Fl (or rather on the 
reduced cross-section at high y) could serve to discriminate more precisely between the 
various alternatives. 

5.2. Varying a s {m 2 ) 

In order to explore the sensitivity of these results to a 8 , we first repeat the analysis 
with a s at the extremes of a range which is representative of the current uncertainty [40] : 
a s (m 2 ) = 0.116 and a s (m 2 ) = 0.122. The x 2 values are shown in Figs. 4a,b while the 
corresponding values of X g are in Figs. 5a,b. For low values of a s both resummations 
give significantly better fits than two loops, while for high values all the best fits are of 
similar quality. We already see at this stage that at two loops the best \ 2 is obtained 
when a s (m 2 ) = 0.122, while in both S and R resummations good fits may be obtained for 
each of the chosen values of a s . The qualitative behaviour of the gluon distribution as a 
function of A is relatively insensitive to the value of a s ; however, as a s increases, all gluon 
distributions become more valencelike to compensate. 

We now describe the results of a second round of fits where also the value of a s (m 2 ) 
is left free. The number of degrees of freedom is now ridof = 92. In fig. 6 we show the 
resulting x 2 as a function of A. Some interesting new features are evident. For most values 
of A, the minimum value of the x 2 generally decreases when compared to the fits with fixed 
a s because we have one more parameter to adjust. The result in the R-resummation is 
similar to that shown in the fixed a s fits, but the range of A where the quality of the fit 
is comparable with that of the unresummed two loop one is now somewhat wider. In the 
S-resummation, the best-fit \ 2 is now always close to the two loop value throughout the 
range of A £ 0.2. Thus, iia s {m 2 z ) is left free parameter, the x 2 is no longer improved 
by implementing the small x corrections. However, if we look at the fitted value of a s (m 2 ), 
plotted in fig. 7, we see that even when the quality of the resummed and unresummed fits 
are essentially the same, the resummed best-fit value of a s (m 2 ) is generally lower than the 
two-loop one. Furthermore, as a function of A, we see that a s {m 2 ) decreases as A increases. 
Indeed, the value a s (m 2 ) ~ 0.122 for the NLO fit is somewhat large in comparison with 
the world average central value a s (m 2 ) ~ 0.119 (as it was in previous fits to older small 
x HERA data [43]). However the resummation tends to bring this value down: in the 
S-resummation in the range A ^ 0, ot s (m 2 ) drops from 0.122 down to a value that can be 
as low as 0.115. For larger values of A the value of a s decreases further but the quality of 
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Figure 4: Same as Fig. 1, but with a s {m 2 z ) = 0.116 (above) or a s (m 2 z ) = 0.122 
(below). 
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Figure 5: Same as Fig. 2, but with a s (m 2 z ) = 0.116 (above) or a s {m 2 z ) = 0.122 
(below). 
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fit deteriorates. This is sufficient to explain why if a s is fixed to the world average central 
value of 0.119, the resummed fits are in better agreement with the data than the two loop 
fit (Fig. 1). In the R-resummation we find again a rapid variation, though a result for 
a s {m 2 ) similar to the two loop one is reproduced in the small range of A where the quality 
of the fit is good. 

In all fits discussed so far, the Q 2 behaviour of A was assumed to be of the form 
^(Q 2 ) = &s(Q 2 )c where c is a constant, independent of Q 2 . However, as discussed in 
sec. 3, lacking a more complete understanding of the effects of running on the value of A 
and on its Q 2 dependence, it is necessary to study the sensitivity of our results to this 
assumption. A possible alternative is to assume that, in the domain of the most relevant 
data, the actual A can be replaced by some average value, independent of Q 2 . We consider 
this as a somewhat drastic assumption; however, it seems appropriate to choose a rather 
extreme alternative given our ignorance of the true behaviour. The fixed A case is also 
interesting for the theoretical reasons discussed in sec. 4: first, in this case there is no 
need to perform an extra subtraction in order to remove the instability of the effective 
anomalous dimension eq. (4.18), and also, the spurious poles (4.22) that appear in the 
effective NLLx BFKL function are absent in this case. We have argued in sec. 4 that the 
effect of both of these problems should in practice be small even when A does depend on 
Q 2 , provided the value of A is itself reasonably small. We can now verify the correctness 
of this conclusion by checking that the results of the fits do not change very much if A 
is taken to be scale-independent. In fig. 8 and 9 we compare the curves of fig. 6 and 7 
with the corresponding ones for A independent of Q 2 . We see that indeed the changes are 
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Figure 7: The best-fit values of a s {m 2 z ) for the fits of Fig. 6. 

indeed quite small, particularly in the case of the S-resummation, thereby supporting the 
reliability of our resummation procedure. 

In fig. 8 and 9 we also show an additional horizontal line: this line corresponds to yet 
another determination of A(Q 2 ), in which we have taken the NLO perturbative result for 
A(Q 2 ), i.e. 

A(Q 2 )=a s (Q 2 ) X o(l/2) + a 2 (Q 2 )xi(l/2). (5.2) 

We see that with this choice, in the S-resummation, the x 2 is essentially the same as the 
two loop one. This is because the perturbative X(Q 2 ) is very small or negative in most of 
the region of the data and this is what is required to fit the data with the S-resummation. 
It is however remarkable that the precise form of the perturbative X(Q) is in such a good 
agreement with the data. Note however that if the R-resummation is adopted instead, the 
perturbative X(Q) is very far from fitting the data because it is much lower than the values 
of A where the fit is acceptable. The dependence of the fitted values of a s {m 2 z ) (Fig. 9) on 
the choice of these different prescriptions for the Q 2 dependence of A is similarly small. 

In order to assess the impact of higher order corrections it is useful to study the renor- 
malization scale dependence of physical observables, such as a s itself. The renormalization 
scale variation is done by letting 



q s (Q 2 ) 
2tt 



<*s(k R Q 2 ) 
2tt 



i + A) 



a s (k R Q 2 ) 
4ir 



Ink 



R 



(5.3) 



everywhere in our computation and suppressing all subleading terms. In the resummed 
computations we take A proportional to a s (knQ 2 ) with the constant of proportionality 
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Figure 9: The best-fit values of a s (m 2 z ) for the fits of Fig. 8. 
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refitted as a function of k^. The ensuing % 2 and a s (m 2 z ) are plotted as a function of kn 
in Fig. 10a, b for the two loop, R-resummation and S-resummation fits. In the two loop 
case, the \ 2 grows quickly when kn is moved away from unity, and a s grows monotonically 
with a particularly rapid rise at large scales. This suggests a large scale uncertainty in the 
value of a s extracted at two loops from HERA data [43,38], due to important higher order 
terms in the perturbation series. By contrast, both the resummed fits remain good even 
if the renormalization scale becomes very large, and the corresponding values of a s are 
surprisingly stable. This suggests that the error due to renormalization scale variation in 
a resummed determination of a s at HERA would be significantly less than in a two loop 
estimate. However, this would have to be offset against the extra uncertainty due to A 
discussed above, as is clear from Fig. 10b. 

6. Summary and Conclusion 

Summarizing, we have compared our small x resummation formalism with structure 
function data from HERA. We recall that our main starting point is to assume that there 
is a region of x and Q 2 , including the HERA kinematic region, where both the leading 
twist Q 2 evolution and the BFKL small x evolution equations are simultaneously valid. 
Using the information from both equations, we construct a double leading expansion for 
the anomalous dimension. We then reorganise the expansion by factorizing out the small 
x behaviour of the splitting function, parametrised in terms of an exponent A. 

A priori, a variety of different resummations are conceivable. Here we discussed two 
distinct possibilities, the S and R resummations, but intermediate solutions could also be 
constructed. We find that the data are in good agreement with the resummed structure 
functions in the S-resummation approach, for a wide range of values of A, starting from 
small positive values down to negative ones. Similarly, in the R-resummation there is also 
good agreement, but here only in a narrow range of positive values of A. In each case 
the range of values of A preferred by the data is such that the anomalous dimension is 
close to the two loop one. This may be evidence that a correct all-order treatment of 
the running coupling in the Regge limit would lead to an anomalous dimension with a 
singularity structure similar to the unresummed two-loop one. 

We think that the data provide a very clear indication that the true value of A is 
small or negative. It is an open question whether or not A can be reliably computed 
in perturbation theory. In this respect it is interesting that in the S-resummation the 
perturbative evaluation of X(a s ) at 2 loops gives a particularly good fit to the data with a 
value of a s {m 2 ) in perfect agreement with the world average. 

An important consequence of our analysis is that even though the NLO QCD for- 
malism leads to an almost unequalled agreement with the data, still the extracted gluon 
density and the value of the strong coupling are significantly affected by the resummation, 
in a way which depends on the precise value of A. 
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n 


a n 


n 


n 


1 


1 





6.191566 


2 





-2.249876 


-3.338781 


3 





0.527795 


-0.551911 


4 


0.112797 


0.224766 


2.271531 


5 





-0.858826 


-1.634828 


6 


0.012658 


0.360812 


0.058893 


7 


0.038170 


0.016245 


1.022436 


8 


0.001601 


-0.417144 


-0.980697 


9 


0.011422 


0.245856 


0.249555 


10 


0.017429 


-0.064702 


0.445905 


11 


0.002607 


-0.206370 


-0.590416 


12 


0.008884 


0.157832 


0.268752 


13 


0.009427 


-0.084159 


0.163843 


14 


0.003002 


-0.097305 


-0.341825 


15 


0.006709 


0.094038 


0.225093 


16 


0.005813 


-0.077607 


0.030933 


17 


0.003016 


-0.041283 


-0.185691 


18 


0.005078 


0.050706 


0.166627 


19 


0.004004 


-0.062323 


-0.024490 


20 


0.002832 


-0.014100 


-0.091471 


21 


0.003904 


0.023129 


0.113295 


22 


0.003016 


-0.046270 


-0.041019 


23 


0.002564 


-0.002371 


-0.037742 


24 


0.003069 


0.006793 


0.071643 


25 


0.002425 


-0.032648 


-0.039719 


26 


0.002277 


0.001495 


-0.009422 


27 


0.002475 


-0.002039 


0.042104 


28 


0.002039 


-0.022298 


-0.031744 


29 


0.002006 


0.001740 


0.003773 


30 


0.002047 


-0.006178 


0.022695 



Table 1: The first thirty coefficients a n , 6°, from which 7 S and 7 SS may be computed. 
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n 




dk 

n 




1 


0.781460 


-0.120225 


0.536650 


2 


0.299133 


0.277452 


1.263035 


3 


0.388067 


0.106613 


0.770611 


4 


0.252563 


0.007350 


0.661812 


5 


0.178383 


0.123681 


0.786193 


6 


0.226323 


0.050700 


0.531599 


7 


0.154733 


0.027507 


0.518416 


8 


0.138915 


0.075887 


0.552254 


9 


0.157449 


0.034088 


0.414474 


10 


0.116020 


0.031678 


0.426128 


11 


0.115329 


0.052513 


0.425327 


12 


0.119970 


0.027497 


0.347879 


13 


0.095680 


0.031055 


0.361618 


14 


0.098463 


0.039209 


0.347915 


15 


0.097079 


0.024334 


0.303998 


16 


0.082914 


0.028822 


0.313664 


17 


0.085594 


0.031073 


0.296830 


18 


0.082037 


0.022438 


0.271865 


19 


0.073841 


0.026192 


0.276778 


20 


0.075504 


0.025864 


0.260999 


21 


0.071575 


0.021013 


0.246641 


22 


0.066829 


0.023660 


0.247790 


23 


0.067477 


0.022391 


0.234531 


24 


0.063926 


0.019769 


0.225970 


25 


0.061122 


0.021414 


0.224628 


26 


0.061018 


0.019973 


0.214079 


27 


0.058081 


0.018611 


0.208590 


28 


0.056334 


0.019498 


0.205838 


29 


0.055760 


0.018207 


0.197663 


30 


0.053434 


0.017520 


0.193751 



Table 2: The first thirty coefficients c n , d„, d^ S with which the quark anomalous dimen- 
sions 7^|, 7H, the Fl coefficient functions C£ , C a L ss and F2 MS coefficient functions 
c lss> C 2,ss ma y be computed. 
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